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ABSTRACT 



Aims. We present an asymptotic theory that describes regular frequency spacings of pressure modes in rapidly rotating stars. 
Methods. We use an asymptotic method based on an approximate solution of the pressure wave equation constructed from a stable 
periodic solution of the ray limit. The approximate solution has a Gaussian envelope around the stable ray, and its quantization yields 
the frequency spectrum. 

Results. We construct semi-analytical formulas for regular frequency spacings and mode spatial distributions of a subclass of pressure 
modes in rapidly rotating stars. The results of these formulas are in good agreement with numerical data for oscillations in polytropic 
stellar models. The regular frequency spacings depend explicitly on internal properties of the star, and their computation for different 
rotation rates gives new insights on the evolution of mode frequencies with rotation. 
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1. Introduction 

The field of asteroseismology has now reached its ag e of matu- 
rity w ith the exploit ation of space miss ions CoRoT (Baglin et alj 
2006) and Kepler (iKoch et alJ l2010h that are gathering stel- 
lar light curves with high accuracy. However, there are still 
unresolved issues that hinder the successful pairing of light 
curve frequencies with pulsation modes, which is crucial to 
obtain detailed information on the inner structure of observed 
stars. One of these issues is the rapid rotation of a star around 
its axis, since the exact nature of rotational effects on pulsa- 
tion mod es is not known. In particular, the centrifugal flatten- 
ing (e.g. iMonnier et al.l (I2007I) ) affects the spectrum of pres- 
sure m odes (p-modes) in a complex way dLignieres & Geo rgeot 
2009). This difficulty mainly concerns non-evolved massive and 
inter mediate-mas s pulsating stars which are typically rapid rota- 
tors ([Rover 2009). Recently though, hints of regular frequency 
spacings have been found in the s pectrum of rapidly rotating 6 
Scuti stars observed wit h CoRoT dGarcia Hernandez eTal . 2009; 
Mantegazza et al. 2012), and this could ease future mode identi- 
fication. 

The recent development of accurate numerical models has 
enabled progress in the comprehension of pulsation modes 
in rapidly rota t ing s t ars. It has been f ound in particular 
dLignieres et al l 120061; iReese et al l 120081 120091) that in the 
rapidly rotating regime a subset of p-modes shows approximate 
regular frequency spacings in the form: 



: A n n + A e { + A m \m\ + a, 



(1) 



where frequencies a> n ^ m are given in the corotating frame. 
Quantum numbers n, I and m correspond to node numbers of 
the mode amplitude distributions, A„, Af and A,„ are frequency 



regularities, and a is a constant term. The approximate for- 
mula in Eq. CD shows a better agreement with numerical re- 
sults towards high-frequencies, thus suggesting that this rela- 
tion is of an asymptotic nature. It should also be noted that, 
from computations of disk-averaging factors, the p-modes fol- 
lowing Eq. CD are expected to be among the most visible ones 
(Lignieres & Georgeot 2009). An example of such a mode can 
be seen in Fig.Q] 

The frequency spacings of Eq. £D are notably simil ar to the 
regul arities described by Tassoul's asymptotic formula dTassoull 
1980) for low degree p-modes in non-rotating stars. Tassoul's 
formula at leading order is 



Mn s + - + - 



with the large frequency separation 



A = 2n 2 



c(r) 



(2) 



(3) 



e-mail: pasek@ir same . ups-tlse . fr 



where c(r) is the radially inhomogeneous sound speed, and R 
the stellar radius. The integer n s is the node number of the radial 
component of the mode, while l s is the degree of the associ- 
ated spherical harmonics, and a s depends on surface properties. 
Tassoul's theory has proved to be very useful for interpre taring 
solar-like oscillations in slowly rotating stars. Indeed, the for- 
mula relates observable quantities, such as the regular frequency 
spacing A, to physical properties of stellar interiors. For rapidly 
rotating stars, it would be clearly desirable to gain insights on the 
underlying physics of the potentially observable regular spacings 
A„, Ac and A m by a similar asymptotic analysis. In this paper 
we derive a formula for these regular frequency spacings in the 
asymptotic regime. 
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Fig. 1. (Colour online) Pressure amplitude P Jd/po on a merid- 
ian plane for a polytropic stellar model, with d the distance to the 
rotation axis and po the equilibrium density. The mode shown 
corresponds to n — 50, £ = 1 and m = 1 at a rotation rate of 
Q/Qk = 0.300, where Q K = (GM/R 1 ^) 1 ' 2 is the limiting rota- 
tion rate for which the centrifugal acceleration equals the grav- 
ity at the equator, M being the stellar mass and R eq the equa- 
torial radius. Colors/grayness denote pressure amplitude, from 
red/gray (maximum positive value) to blue/black (minimum neg- 
ative value) through white (null value). The thick black line is the 
ray y located in the center of the main stable island. 



The generalization of the p-mode asymptotic theory to 
rapidly rotating stars is not trivial. Tassoul's theory requires sep- 
aration of variables, which is no longer possible when the star is 
flattened by rotation. For non-separable wave systems, a well- 
known technique to obtain eigenmodes is to study the short- 
wavelength limit of the propagating waves. This limit gives an 
equation for the propagation of rays that is similar to the geo- 
metrical optics limit of electromagnetism, or the classical limit 
in quantum mechanics. Then, by imposing quantization condi- 
tions on the phase of waves propagating on these rays, one ob- 
tains the eigenmodes of the wave system. This technique was 
first developed in the context of quantum physics, and is often 
called semiclassical quantization. 

For spherical stars, the ray limit of pressure waves has 
been previously used to recover the Tassoul asymptotic formula 
from the Einstein -Brillo uin-Keller (EBK) quantization of ray 
dynamics dGoughll"l993l) . This analytical approach is possible 
only when the ray system is integrable. A dynamical system 
is said to be integrable when it has as many conserved quan- 
tities (energy, angular momentum, etc.) as degrees of freedom 
(IOttll2002h . In rapidly rotating stars, there are not enough con- 
served quantities to ensure inte grability of the ray dynamics. 
Indeed, in Lig nieres & Georgeodl200l 12009), it has been found 
that acoustic rays in rotating stars have a very different dynam- 
ical behavior depending on their initial conditions in position- 
momentum space (the so-called phase space). For a polytropic 
stellar model, the numerical integration of the equations for 
acoustic rays displayed various types of solutions. Indeed, one 
can obtain either stable rays staying on torus-shaped surfaces 
in phase space which form structures such as stable islands, or 
chaotic ray s that are dense and ergodic on a phase space volume 
dOttll2002l) . 

A similar behavior has been found in many systems studied 
in the field of theoretical physics known as quantum chaos or 



wave chaos dGutz willed 1 1 990h . This field has among its objec- 
tives to analyze quantum (resp. wave) systems whose classical 
(resp. short-wavelength) limit is partly or fully chaotic. In this 
framework, one can predict the existence of some eigenfunc- 
tions (resp. mode amplitudes) and energies (resp. frequencies) 
of the quantum (resp. wave) system from the different s tructures 
that are present in the classi cal system phase space dPercivall 
19731: iBerrv & Robnikl ll984l). In t he stellar pulsation setting, 
Lignieres & GeorgeotN2008l l2009h found that the mixed (i.e. 
regular and chaotic) character of the acoustic ray dynamics in 
rapidly rotating stars results in a classification of p-modes in two 
broad families: regular modes either associated with stable is- 
lands or whispering gallery zones, and chaotic modes associ- 
ated with ergodic regions in phase space. For the regular modes 
associated with stable islands, the so-called island modes, it is 
known to be possible to obtain approximate analytical solutions 
by solvin g the wave equation i n the vicinity of a periodic stable 
ray dBabich & Buldvrevlll991l) . A simple application of such a 
method is found in modes of optical resonators, where the peri- 
odic stable light ray is a s traight line between two reflecting mir- 
rors dKogelnik & Li 1966). These methods have bee n previously 
emplo yed to obtain m odes of more complex lasing dTureci et al.l 

120081) cavities as well as 
In this paper, 



2002) and electronic (Zalipa ev et al 



quantum chaos systems ( Vagov et al 



2009). 



we 



apply this approach to rapidly rotating stars. 

In the present analysis, we thus construct an asymptotic for- 
mula for regularities in the p-mode spectrum of rapidly rotating 
stars. Part of th e results were alre ady presented in the short com- 
munication of Pase k et al] d201 lh . In the present paper we give a 
detailed derivation of these results, specify their domain of valid- 
ity, extend them with a study of rotational splittings, and explore 
their astrophysical applications. 

The paper is organized as follows. In Sect. |2]we present the 
wave equation for p-modes in rotating stars and its asymptotic 
limit leading to an equation for acoustic rays. In Sect. [3] we 
use a stable periodic solution of the ray dynamics to obtain a 
semi-analytical formula for the associated p-modes, and to de- 
rive a formula for the associated regular frequency spacing s. We 
then compare the results obtained from the derived formulas for 
mode frequencies and spatial distributions with numerical re- 
sults (Sect. HJ. Finally, we suggest directions on how these re- 
sults could be used for the asteroseismic diagnosis of rapidly 
rotating stars by discussing the phenomenological implications 
of the theory in Sect. [5] 



2. P-modes in rotating stars and their asymptotic 
limit 



In Sect. l2.ll we introduce the wave equation for p-modes in rotat- 
ing stars. We then present the asymptotic limit of this equation 
in order to obtain an equation for the dynamics of acoustic rays 
(Sect. El). 

2.1. Pressure modes in rotating stars 

We start with the equation for small adiabatic time-harmonic 
perturbations of the pressure field in a self-gravitating gas. Since 
we are interested in obtaining an asymptotic theory for p-modes 
in the high-frequency regime, we use the Cowling approxima- 
tion (i.e. we neglect the perturbations of the gravitational po- 
tential), an approximation known to be valid for high- frequency 
perturbations in non-rotating stars dAerts et al.ll2010i ). We also 
neglect the Coriolis force. Indeed, in the high-frequency regime, 
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the time scale associated with this force is much longer than 
the mode period, and thus the influence of the Coriolis force 
on pulsatio n frequencies is weak . This has been numericall y 
checked in iLignieres et all (120061) : iReese et alJ (120061 l2008h . 
In the asymptotic regime of p-modes, the oscillation frequen- 
cies are far greater than the Bmnt-Vaisala frequency and thus 
we can discard the terms corresponding to gravity waves. With 
these assumptions, the equation for pressure perturbations is a 
Helmholtz equation such that 



A*P + 



-T = 0. 



(4) 



where *P = P/f is the complex amplitude associated with the 
pressure perturbation P - Re[P exp(-icot)], f is a function of 
the background model, co c is the cut-off frequency of the model 
and c s its inhomogeneo us sound velocity (for a deta iled deriva- 
tion of this equation see lLignieres & Georgeotl2009l) . The stellar 
model is not spherically symmetric due to the centrifugal dis- 
tortion, but is however cylindrically symmetric with respect to 
the rotation axis. Therefore, we can write the pressure field as 
¥ = ¥ m exp(imip) where m is an integer and <p is the azimuth an- 
gle of spherical coordinates. By inserting this expression in Eq. 
© we obtain (cf. Sect.© 



Ad>„, + 4 



2 2 
to - CO - 



d 2 



(5) 



where d is the distance to the rotation axis. The new mode ampli- 
tude <S> m is such that O,,, = yfd¥ m . We introduce a renormalized 
sound velocity: 



c s = 



(6) 



We notice that besides its spatial dependence, c s also depends 
on at and m, and that m is taken as a parameter for the two- 
dimensional wave equation Eq. (|5). 

2.2. Ray limit of p-modes 

In non-rotating stars, the asymptotic theory of high freq uency 
p-mode s has first been derived by IVandakurovl (119671) . and 
iTassoull (Q.980). The method was to use the spherical symmetry 
of the star to reduce the problem to a one-dimensional equa- 
tion in order to obtain the mode frequencies. This method is 
not applicable when the centrifugal force breaks the spherical 
symmetry of the star. In this case though, one can study the 
short-wavelength limit (u> — » oo) of the wav e equation Eq. 
(01) (as detailed in Ligni eres & Georgeotl d20 09)). This provides 
a Hamiltonian system describing the pr opagation of acoustic 
rays. T he Hamiltonian has been derived in lLignieres & Georgeotl 
(2009|)as 



H = - 



2c? 



2 2 2 \ 

or c (r s mr \ 



co 2 d 2 



(7) 



where the frequency-scaled wavevector k p is the projection of 
k = k/o) onto the meridional plane of the star. We notice that 
this expression has been derived from the short-wavelength limit 
of the three dimensional wave equation Eq. (|4]) and then, pro- 
jected onto the corotating meridian plane. An alternative deriva- 
tion would be to start from the two-dimensional wave equation 
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Fig. 2. (Colour online) Poincare Surface of Section (PSS) at the 
rotation rate Q/Q^ = 0.589 for quantum number m = 0. Each 
dot corresponds to the crossing of an acoustic ray with the equa- 
torial half -plane in the (r/R eq , k r /u>) phase space, r being the ra- 
dial coordinate and k r the associated momentum. R eq is the equa- 
torial radius and cj the mode frequency. Red/dark gray denotes a 
chaotic ray, green/light gray a whispering gallery ray, blue/black 
a stable island ray (see text). Upper inset is a close-up of the 
main stable island. 



Eq. (0. In this case, the ray limit yields the same expression with 
the addition of the 1/4 factor of Eq. (|5j that accounts for the im- 
possibility of acoustic rays to go through the rotation axis (i.e. 
d = 0). Throughout the paper we use Eq. (0 as the Hamiltonian 
for acoustic rays. 

To probe the integrability property of a dynamical system, 
it is convenient to use the Poincare surface of section (PSS) , 
a standar d tool in dynamical systems theory (lGutzwillerl ll990: 
|Ot3|2002) to visualize the structures in phase space. A PSS is a 
lower dimensional slice of phase space. The acoustic ray dynam- 
ical system in the meridional plane has two degrees of freedom, 
which gives a phase space of dimension four (two for positions, 
and two for momenta). There is one conserved quantity in the 
form of the acoustic wave frequency, so the dynamics belongs 
to a three-dimensional manifold in phase space. By fixing an 
additional position or momentum coordinate, we obtain a two- 
dimensional PSS which can be easily visualized. An example 
of such a section for our system is shown in Fig. [2] Different 
choices o f PSS variables are possible , some of which are pre- 
sented in lLignieres & Georgeotl (12009). We have here chosen to 
fix the colatitude 9 = jt/2, so that the PSS corresponds to the 
crossing of rays with the equatorial half-plane. We thus display 
a section in coordinates (r/R eq , k r /u>) where k r is the norm of the 
radial wavevector, o> the mode frequency, and R eq the equatorial 
radius (that may be greater than the polar radius since the star is 
flattened by rotation). In such plots, each dot corresponds to the 
crossing of an acoustic ray with the PSS. Successive dots from 
a single ray will form lines in integrable zones, or fill surfaces 
densely in chaotic zones. We see in Fig.|2]that when the rotation 
rate Q./Q.K (where Q. K = (GM/R 3 eq ) 1/2 is the limiting rotation 
rate) is large, different structures coexist in the system phase 
space: stable islands correspond to concentric circles around 
a stable ray, whispering gallery rays to lines near the surface, 
and chaotic zones to densel y filled areas (for more details, see 
ILignieres & Georgeotll2009l) . In this paper, we will focus on the 
2-periodic stable island which is the main stable island (shown 
in the inset of Fig.|2]i. The rays' dynamics is very sensitive to the 
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rotation rate, so the PSS will be different for each rotational ve- 
locity. Indeed, the locus in phase space of the main stable island 
changes as rotation increases. The major change happens when 
the central ray of the 2-periodic stable island undergoes a bifur- 
cation at Q./Q.K — 0.26. For slow rotation rates, the central ray of 
the island is located on the polar axis, and through this bifurca- 
tion it transforms into two stable rays surrounding one unstable 
ray on the polar axis. Then, as rotation increases, the stable is- 
land will coast away from the polar axis. This bifurcation will be 
of some importance in the following, when we will show that one 
can construct approximate eigenmodes of the wave system from 
this 2-periodic stable island. In Fig. [1] one can see an example 
of an island mode obtained from a full-numerical computation, 
together with the central ray of the 2-periodic stable island for 
the same rotation rate and quantum number m. 



3. Semi-analytical method for island modes 

In this section we construct an asymptotic approximation of a 
subset of regular p-modes associated with a stable periodic ray. 
The method is based on th e works of Babich and coworkers (see 
iBabich & Buldvre"vl(|1991l) and refer ences therein for the general 
formalism, and e.g. IZalipaev et al.l d2008l) : IVagov et"aD {2009) 
for some applications). It consists in deriving an approximation 
of the wave equation in the vicinity of the ray (Sect. [XT), finding 
Gaussian wavepacket solutions related to the stability properties 
of the ray (Sect. I3.2K and then deriving the asymptotic formula 
for the frequencies from a quantization condition (Sect. I3.3l l. 



3.1. Approximate wave equation in the vicinity of a stable ray 

For a given rotation rate and quantum number m, we start with 
the central periodic ray of the main stable island (see Sect. 
13. This ray must be computed by numerically evaluating the 
Hamiltonian equations derived from Eq. ((7}. In the following we 
will call this ray y. The first step is to write the wave equation 
Eq. (0 in the vicinity of y. For this, we use a local orthonormal 
coordinate system (s, defined as r — sT+^N where s is the arc 
length along the ray, T the unit tangent vector, £ the transverse 
coordinate and N the unit vector normal to T. The two basis vec- 
tors are related by the curvature k(s) of the ray as follows: 



1 



<s) = -— ■ T . 

ds 

In this coordinate system, the wave equation Eq. (0 reads 



(8) 



l Id lhd<b m \ d (h s d<s> m \\ u> 2 
h s ht \ds\h s ds J + % [hf d% J) + c s (s,o 2 "' ' 

where the scale factors h s and h s (lArfken & Weberll2005l) are 

h 2 s =(T-tK(s)T) 2 = (l-f<c(s)) 2 (10) 



and 



hi = N(s) 2 = 1 



(ID 



In the vicinity of the ray y, that is for small the terms of Eq. 
(0 are simplified using: 



1 



l+^( S )+fK(s) 2 + 0(f), (12) 



(W*C0) 
1 



1 + 2%k(s) + 3fK(s) 2 + 0(f) , (13) 



1 + 3£k(s) + bf-Kisf + 0(f) , 



(14) 



and 



1 



1 



c s (s,f) 2 c s (s,0) 2 



lf=0 



{ + 0(f). (15) 



We then express the function O m (s, in terms of a WKB ansatz 
as 

O m (s, £) = expOW)£/,„(s, £ co) , (16) 

where t is an unknown function of position. The fundamental 
assumption underlying the theory of Babich is that, asw-> +oo, 
the mode is localized on the acoustic ray and that its transverse 
extent scales as 1 / -\fa>. Such a solution can be found by assum- 
ing that the transverse variable £ scales as 



f = o(i/V5) 



(17) 



Then, from an expansion of Eq. (0 in u>, one obtains at the dom- 
inant order that the WKB phase in Eq. ( fT6b depends only on s as 
dr - ds/c s . At the next order in at one finds a parabolic equation 
for the function V m : 



SrVm 
dv 2 



- K(s)v 2 V m + 



2i dV m 
c s (s) ds 



= 0. 



with 



K(s) 



1 d 2 c, 



c s (sy df k=o 



(18) 



(19) 



where we introduced the scaled coordinate v = £ and V,„ = 

Uml a/Q- 

3.2. Solutions of the parabolic wave equation 

To find a solution to Eq. ( fT8l , we first find a solution at a fixed 
arc length s, and then study how this solution must evolve with 
s. At fixed s, the first two terms of Eq. (fT8l correspond to the 
equation for a quantum harmonic oscillator in the direction N 
transverse to y. Thus, as we know from quantum mechanics 
dCohen-Tannoudii et alj[l973l) . a solution of this equation is a 
Gaussian wavepacket, transverse to the ray, that we write as 



v:: = ,i(.v)cx P |i^ 



(20) 



with r an unknown complex-valued function. To find the varia- 
tion of this Gaussian wavepacket along the ray y, we introduce 
a solution of this form in the parabolic equation Eq. (fT8~t and 
obtain a Riccati equation for T 



c, ds 



and a simple form for the factor A 
1 dA 

A ds ~ 2 



(21) 



(22) 



In the following, we show that the equation for F is related to the 
ray properties in the vicinity of y and can thus be solved from 
ray dynamics computations. First, using the variables (z(s), p(s)) 
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defined as T(s) = ± and p{s) = l^, Eq. <ED is trans- 
formed into the Hamiltonian system: 



dz 

dr 
dp 

dr 



-c 2 ,Kz ■ 



where the (time-dependent) Hamiltonian function is 



2 2 

H (p,z,r) = c 2 s ^-+c 2 s K-, 



(23) 
(24) 

(25) 



t being the time coordinate. Equations ( 1231241 have two in- 
dependent solutions (z, p) and (z, p) that are conjugate to each 
other. Of these two solutions, only one is physically relevant, i.e. 
corresponds to a localized wavepacket. According to Eq. (l20l > 
this happens if Im(F) > for all s, and as shown in Sect. [B] the 
sign of Im(F) stays constant along y since 



Im(T) 



1 1 

2W 



(26) 



The variation along y of the Gaussian wavepacket can now be 
linked to the dynamics of the acoustic rays nearby y. Indeed, 
Im(r) has a simple expression in terms of the complex variable 
z, as shown in Eq. (f26b . This variable, on the other hand, obeys 
Eqs. d23H24t which can be shown to be the same as the equation 
describing the deviation from y of a ray nearby y (see deriva- 
tion in Sect. Q. The Hamiltonian in Eq. d25b is thus a local in- 
tegrable approximatio n, also known as a normal form approxi- 
mation dAmpl'dll 19891) . to the full Hamiltonian for acoustic rays 
written in Eq. (|7). 

Now, our task is to find the two linearly independent complex 
conjugate solutions of Eqs. ( I23ll24b . The terms in these equations 
depend only on quantities that are evaluated on the periodic ray 
y. Therefore, these equations are periodic in s, or equivalently in 
r. Eqs. ( l23ll2"4T > can thus be written as: 



where the matrix S 



d_lz 
dr \P 



S(T) 



S(T) 



c 2 
-c 2 K 



(27) 



(28) 



verifies E(t + T y ) = S(r), T y — <j5 if being the time period as- 
sociated with y. Then if (z(t), p(r)) is a solution of Eq. ( l27l i. so 
is (z(t + T y ), p(r + T y )) and the two solutions are related by the 
following linear map 



Z(T + Ty) 
P(T + Ty) 



M 



z(t) 
P(r) 



(29) 



where M is called the monodromy matrix (ICvitanovic et al.1 
(1201 Ol) . and references therein). As (z, p) describe ray devia- 
tions from y, the matrix M characterizes the stability of y. As 
y is stable, we know that |Tr(M)| < 2 and that the eigenvalues 
are of modulus one and complex conjugates of each other i.e. 
A* = exp(±/a') with a e]0, n[ (cf. Sect. iDt . where a is called a 
Floquet phase or stability angle. Hence the two linearly indepen- 
dent solutions of Eq. ( 1271 ) can be written in the form: 



(z(t),p(t)) = exp|±z'— t|k ± (t)v 



(30) 



are 



where the functions u±(r) are periodic with period T y and v ± 
independent eigenvectors of the monodromy matrix M. 

An expression of the monodromy matrix in ter ms of second 
deriva tives of the action function S can be derived dBogomolnvl 
2006). The action function S is defined by a trajector y from the 
positi on qi to qf for a given energy or frequency to dGutzwillerl 
Il990h . For our purposes, the action is written as 



rif 1 

S(qi,qf,cS)- I —dcr, 

J a, c s 



(31) 



where cr is the arclength along a ray nearby y. If we write the 
monodromy matrix as 



Zf\ = [Mn M n \lz, 
Pf] \M 2 \ M 22 



(32) 



we can express its components from the second derivatives of 
the action function S with the following formulas 



d 2 s 

dzidzj 



1 



dz 1 , 



M, 2 



d 2 S 



dz} 



M 12 
M, 2 



(33) 



where the positions qi and qf are written as (suzd an d (• s /»Z/)> Zi 
and Zf being respectively the initial and final transverse positions 
of the neighboring ray after one period, and the derivatives are 
evaluated on the periodic ray. From these expressions, and the 
simple formula giving the roots of a second degree polynomial 
(cf. Sect. [D]), we can obtain an expression for the stability angle 
a as 



a = arctan 



' V-Tr(M) 2 +4' 
Tr(M) 



where 



Tr(M) = - 1 



/ d 2 s r 1 


'd 2 s 


d 2 s' 


\dzidzf) 


M 


+ dz 2 



(34) 



(35) 



We thus have obtained a solution of the approximate wave 
equation Eq. ( fT~8b in the form of a Gaussian wavepacket (Eq. 
d20li). whose evolution along the ray y is given by Eqs. d27ll28l i. 
It is possible to find other solutions of Eq. dT8l that have a finite 
number of nodes in the directi on transverse to y. In the sam e way 
as in quantum mechanics dCohen-Tannoudii et al.|[l973l) . these 
solutions can be obtained from Eq. d20T l using the annihilation a 
and creation operator a 1 that are defined as 



• 9 H-f - 3 ~ 

a — —iz— pv and a = — iz- pv . 



(36) 



3v dv 
These operators have the commutation rule [a, a* J = 1 , where 

the commutator is defined as [A, fij = AB - BA. It can be shown 
that the expression for the higher-order solutions is 



Vi 



(37) 



This defines a recurrence relation, whose solutions are propor- 
tional to Hermite-Gauss polynomials and can be written in the 
following way 



^-(j) w ft(VEa>)=g3 



(38) 



with He the Hermite polynomials of order I. Finally, we can 
write the solutions of Eq. (0 as 

<(5,v)= y[F s V e m (s,v)cx V (ioj j . (39) 
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3.3. Quantization condition and regular frequency spacings 

The quantization condition is based on the single-valuedness of 
the solution presented in Eq. d39l >. and thus asserts that the phase 
accumulated by the function <J>£ over one period must be a mul- 
tiple of 2tt. In the following, we assume without loss of gener- 
ality that the eigenvalue which corresponds to the wavepacket 
localization is exp(+/aO. The contribution of function V e m to the 
dynamical phase of <E>£ is obtained from Eq. (l30t and Eq. ( l38l l. 
We obtain that the phase accumulated over one period is 



X 1 , ° 

■>n,l,m (D — ds 

Jy Cs 



+ 2nN r 



- (a + 2nN r )£ = 2nn + x , (40) 



where T y = ^ jr is the acoustic travel time along y and a the 
Floquet phase that is defined modulo 2n. For our purposes, we 
must also take into account the number N r of multiples of 2n 
acquired by the phase. N r can be computed by following the 
evolution of the eigenvector v over one period, and we verified 
numerically that, alternatively, A^,- is also the winding number 
around y of a ray nearby y d uring one period . The last term in 
Eq. ( f4Qb is the Maslov phase dGutz willerll 199 0) that comes from 
the reflection of the wave on the boundaries. The formula for the 
frequencies, co„j, m , of island modes is thus 



1 



lis 



2/r|n+M + k+i|(2/r ; V, f <r) 



or, in a form that makes the regularities more visible, 

w„,/,„, = 6 n (m)n + 6 e {m)€ + j3(m) , 



with the frequency regularities 
2n 



6„ = 



and the constant term 



and 6( = 



2nN r 



Jy C s 



5 n + 5[ 



(41) 



(42) 



(43) 



(44) 



The quantum numbers n and t correspond to node numbers of 
the p-mode, respectively in the longitudinal and transverse di- 
rections of the central ray y as illustrated in Fig. [3] This is to 
be contrasted with the case of a spherical mode where the most 
natural labeling are the quantum numbers of spherical harmon- 
ics. Eq. d4TT > is a semi-analytical formula since the quantities T y , 
N r and a must be computed numerically from the Runge-Kutta 
integration of the Hamiltonian equations for acoustic rays. The 
acoustic time, T y , is directly computed from y itself. From the 
intersections of a ray nearby y with the PSS (as can be seen in 
Fig.0, we compute the monodromy matrix M that maps one in- 
tersection with the PSS to the next one. Then, by diagonalizing 
this matrix one obtains the Floquet phase a from its eigenvalues, 
and the functions z and F from its eigenvectors. It is important to 
note that for m even, only modes symmetric with respect to the 
rotation axis exist. Since the preceding theory does not take this 
phenomenon into account, the theoretical value of 5( is multi- 
plied by two when the ray y coincides with the rotation axis, i.e. 
for rotation rates less than the bifurcation point. Finally, it can 
be noted that a formula similar to Eq. (HTt can also be obtained 
through the fo rmalism of th e Gutzwiller trace formula following 
the method of Millet] dl975l) . 




m = 3 
m = 2 
m=l 
in < 
m = - 1 
m = -2 
m = -3 



ii = 2n 

e=o 

n = 2n +1 



n = 2n +1 



| Even 
I Odd 



I = 



= 1 I =2 1=3 



Fig. 3. (Colour online) Left: Schematic representation of the 
mode labeling/quantum numbers used for island modes. Right: 
Illustration of the relation between the quantum numbers of 
spherical and island modes. The quantum number n is the num- 
ber of nodes in the longitudinal direction of the island mode 
(along the ray y), I is the number of nodes in the transverse di- 
rection of the mode (transverse to y) and m is the number of az- 
imuthal nodes. n s , £ s and m s are the quantum numbers of spher- 
ical modes. In the right figure, the multiplets of island modes 
correspond to diagonal colored bands, whereas the multiplets of 
spherical modes would have the form of vertical bands. The plot- 
ted mode on the left corresponds to n — 46, 1 — 1 and m = 0. 



4. Comparison with numerical results 

As the present asymptotic theory relies on various assumptions, 
its relevance for stellar seismology is not guaranteed and needs 
to be assessed through a comparison with exact calculations of 
realistic stellar models. In this section, the comparison is done 
with numerically computed modes in uniformly rotating poly- 
tropic models of stars. The hypotheses of the asymptotic theory 
are the following : first, it is valid in the asymptotic regime, that 
is, for high enough frequencies. Second, the island modes are 
constructed from a stable island of acoustic ray phase space. At 
null rotation such a structure does not exist, so we expect that the 
theory fails to describe spherical mode amplitudes and frequen- 
cies. A stable island immediately appears at non-zero rotation, 
but its phase space volume must be high enough for an island 
mode to exist. This volume increases with f requency and rota- 
tion Csee lLignieres & Georgeoj (12008. 2009) for details). Thus, 
for a given frequency range, the number of island modes starts 
from zero at small rotation rates and progressively increases as 
the rotation and thus the phase space volume of the stable island 
grows. Actually, low degree spherical modes become progres- 
sively island modes as rotation increases. Another assumption 
used in finding a solution to the wave equation is that the mode 
decays as oc 1/ -\Ja> in the direction transverse to the periodic 
ray. Finally, the theory also neglects the Coriolis force and the 
perturbations of the gravitational potential. 

In the following, the asymptotic theory is compared with 
highly accurate computations of high frequency adiabatic modes 
in uniformly rotating polytropic stellar models with index N = 3, 
the Coriolis effect and perturbations of the gravitational potential 
being taken into ac count. The accurac y of these calculations, de- 
scribed in detail in Rees e et al . (2006), is very high (the relative 
precision on the frequencies is 10~ 7 ) and thus does not interfere 
with the present comparison. A large number of modes were fol- 
lowed from Q./Q.K — to Q/Q^ = 0.896. At zero rotation these 
modes are low degree ( s e {0, 1,2,3), high order n s e [21,25] 
modes. At higher rotation rates, they become island modes and 
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Table 1. Island mode quantum numbers n, (, m of numerically 
computed modes, corresponding to n s e [21,25],£ s e {0, 1,2,3), 
m s e [—6g, € s ] in terms of spherical mode quantum numbers. 



n 


I 


m 


42, 44, 46, 48, 50 





-3, -2, . . . , 3 


42, 44, 46, 48, 50 


1 


-1,0, 1 


43, 45, 47, 49, 51 





-2, -1 2 


43, 45,47, 49,51 


1 






can thus be labeled with n and I, the number of nodes along and 
transverse to y, respectively, as illustrated in Fig. [3] The relation 
between the quantum numbers at zero and high rotation rates is 
dReesell2008h : 



n — 2n s + [(( s + m s ) mod 2] , 

_ €„ - \m,\ - [(( s + m s ) mod 2] 
€ ~ 2 
m — m* . 



(45) 
(46) 
(47) 



We remind here the reader that rotational multiplets are defined, 
in the non-rotating case, as a set of frequencies with identical n s , 
( s quantum numbers but different values of m s for m s e [-f s , (A- 
For rotating stars, we can define multiplets as frequencies with 
identical n and i but different m e Z, i.e. without any limiting 
value for m. The relation between the two sets of quantum num- 
bers and the two types of multiplets is illustrated in Fig. [3] In this 
figure the multiplets of island modes correspond to diagonal col- 
ored bands, whereas the multiplets at zero rotation would have 
the form of vertical bands. We restricted ourselves to numeri- 
cal modes with ^™ ax = 3 so, in terms of island mode quantum 
numbers, the range of numerically computed modes is the one 
given in Table Q] The associated numerical frequency spacings 
are defined as: 

(48) 



' J n+\,l,m ' 



°n,e,m 



and 



J n,C+\,m 



(49) 



The semi-analytical asymptotic theory also requires determin- 
ing the a term in Eq. d43l numerically. To test the robustness of 
this calculation, we checked that the frequency spacing 6( only 
weakly depends on the choice of the ray nearby y that is used 
to compute a. Also, the spacings 6„ and S{ neither depend on 
the resolution of the background model nor on the integration 
parameters of the Runge-Kutta method. 




0.2 0.3 0.4 
Q./Q.K 



0.6 



1.2 
1.0 
0.8 
(Up 0.6 
0.4 
0.2 
0.0 



A 



A 

. 1 • • ■ i • • 



i 



A A A A A A 



0.0 



0.1 



0.2 



0.3 
Q./Q.K 



0.4 



0.5 



0.6 



Fig. 4. (Colour online) Comparison of frequency spacings S n j 
between island modes, computed from numerical simulations 
and semi-analytical formulas for different values of Q./Q.K (the 

frequency spacings are normalized by u> p with 

R p the polar radius). Circles: 6„, triangles: S(, red/dark gray: 
semi-analytical results, blue/black: numerical results. Numerical 
results correspond to different sets of S„ and 5^ values, with 
n e [42,51], t e {0, 1} for m = 0; and n e [42,51], t = 0, 
n s {42,44,46,48,50}, i = 1 for m e {-1,1}. Upper panel: 
m = 0. Lower panel: m = ±1. 



4. 1 . Regular frequency spacings 

According to Eq. (l42l . the structure of the spectrum is character- 
ized by the two spacings S n (m) and S((m). In Fig.|4j their semi- 
analytical and numerical values computed for m = and \m\ = 1 
are compared as a function of the rotation rate. One can see that 
the semi-analytical regularities 6„, 5f, and the full computations 
of high-frequency p-modes are in good agreement for almost all 
rotation rates. For m = 0, around Q./Q.K 0.26, the agreement 
degrades significantly. In this rotation range, the ray y in the cen- 
ter of the main stable island undergoes a bifurcation from one 
stable ray on the polar axis to two stable rays surrounding one 
unstable ray. When such a bifurcation occurs, the eigenvalues 
of the monodromy matrix become A* = 1, corresp onding to a 
Floquet phase a = mod 2n (see e.g. Brack 2001). This is in- 
deed what happens at Q./Q.K — 0.26, as S( °c a goes to zero. 
Such a behavior conveys the non-validity of the present normal 



form approximation for rays undergoing a bifurcation. One pos- 
sibility would be to use other local appro ximations of the ray 
dynam ics called uniform approximations (ISchomerus & Sieberl 
1997). The discrepancy coming from the bifurcation is not to be 
found for m + 0, since the stable ray stays away from the polar 
axis and does not undergo a bifurcation as rotation increases. 

Although, as mentioned before, the theoretical and numeri- 
cal frequency spacings are not expected to match for slow rota- 
tion rates, the discrepancies remain small in this rotation range. 
This is due to the fact that, as O/Qjf approaches zero, the sta- 
ble ray is along the polar axis and, according to the expression 
of S„, this implies that S„ - A/2, i.e. half the large separation 
defined in Eq. (0. Now, using the first order of Tassoul's for- 
mula and the quantum numbers conversion rules Eqs. (I45H47K 
it is easy to see that, at zero rotation, S„ = a> N , , , - w", 
is expected to be close to half the large separation. Note also 
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Fig. 5. (Colour online) Comparison of frequency spacing 6 m be- 
tween island modes, computed from numerical simulations and 
semi-analytical formula Eq. ( |5TT i for different values of Q/£2# 

(the frequency spacings are normalized by to p with 

R p the polar radius). Red/dark gray: semi-analytical results of 
Eq. ( |5TT i for n — 42 and m — 1, blue/black: numerical re- 
sults. The numerical results correspond to different sets of 6 m 
forne [42,51], £ = 0andm e{-l,l). 



that the doubling of the numerical values observed in Fig. |4] at 
small rotation rates is due to the small separation that appears 
at the next order of Tassoul's theory. Concerning 6c, the semi- 
analytical calculations indicates that 6c goes to 26„, that is A, for 
slow rotation rates. Again, Tassoul's formula applied to the n, 
i quantum numbers shows that 6 N t is close to A. For these rea- 
sons, the frequency spacings 6„ and 6c converge to the results 
of the first order of Tassoul's formula, though their derivation is 
formally not possible for non-rotating stars. 

In order to investigate the drift between the spectra of differ- 
ent m as rotation increases, we consider the frequency spacing: 



(50) 



Figure|5]displays a comparison between the numerical and semi- 
analytical values of 6„, for I — and |m] = 1. As expected, the 
agreement is not good at small rotation rates. Using the first or- 
der of Tassoul's formula to approximate the numerical results at 
zero rotation, 6% = co N , - to N is found to be close to |ot|A/2 
when Q./Q.K = 0. This is not compatible with the asymptotic the- 
ory of the island mode since it predicts that 6 m goes to zero when 
to goes to infinity. Indeed, 6„{m) depends on m/co because c s and 
the ray path y, given by the Hamiltonian Eq. (Q), both depend 
on m/co. An alternative explanation is to consider the spatial dis- 
tribution of island modes of fixed m and (: one finds that in- 
creasing to produces both larger derivatives along the stable ray 
associated with a higher node number n, and larger transverse 
derivatives because the transverse extent scales as 1 / ^[to. Thus, 
the contribution of the azimuthal derivatives becomes negligible 
in the wave equation Eq. ©. We also verified that 6,„ displayed 
in Fig.0diminishes when n is increased. Thus, for rotation rates 
such that the numerical modes are not fully island modes, they 
behave more like spherical modes and 6% shows clear discrep- 
ancies with the asymptotic results. By contrast, at high rotation 
rates, an approximate analytical formula for 6 m is derived in the 
following and shown to closely reproduce the numerical results. 



Starting from 

6 m ({ = 0) = [6„{m) - 6 n (0)]n + \J3(m) - p(Q)) , (51) 

we first assume that n is large enough to neglect f3(m) - 0(0). 
From Eq. (l43t . 6„(m) is equal to 2nlT y {m) where 



l ~~ 2 



d 2 



(52) 



The dependence of T y in m is explicit in the integrand but is 
implicit in the integration path y. In the following, the variation 
of the location of y with m/to is assumed to be negligible. Then, 
an expansion in l/a> of the integrand in Eq. d52l leads to: 



c s (m 2 - 1/4) 
J 2 



ds 



(53) 



Hence we obtain an approximate expression for 6„(m) of the 
form: 



6„(m) =a 



27T It §y f dS + §y 



c.,(m 2 -l/4) 



d 2 



ds 



Jy C s 



(54) 



If we insert the previous expression for 6 n (m) in Eq. ( BTT l and 
neglect f3(m) - yS(0), we have 



6 m {( = 0) - 



(4?) 



(55) 



Finally, normalizing by to p = ^GM/Rp and replacing wbyn^ 
yields 



S m {t = 0) 



(m_\ 2 $% ds I n \ 2 J_ 



I m 



(56) 



\y/n) 4naj p J y d 2 



where we have used the fact that — j— stays constant in the 

a/Up J 

frequency range considered here, and is known to be close to 
cop ^-/2n. The previous expression will be made more precise 

by renormalizing the value of c s by y/l — a> 2 /to 2 to take into 
account that co c is not negligible, and indeed is of the order of 
to, close to the stellar surface. In Fig. [6] the numerical values of 
5 m (t — 0) as well as results for the last term in Eq. (l56l are plot- 
ted as a function of ml V« for £2/Q# = 0.419, showing a good 
agreement. This behavior is valid for rotation rates higher than 
Q./Q.K — 0.4. It must also b e noted that in the numerical calcula- 
tions bv lReese et al. (2009), using more realistic stellar models, 
the asymptotic -^j= dependency was also found empirically. 



4.2. Pressure amplitudes of island modes 

In this section, we compare the results obtained from the semi- 
analytical formula for mode spatial distributions Eq. ((39) with 
results from full numerical computations. Equatorial cuts of 
the semi-analytical modes can be expressed as a function of 
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Fig. 7. (Colour online) Examples of normalized amplitudes distributions (real part of 0£) on the equator as a function of position 
r/R e „ (where R eq is the equatorial radius of the stellar model) for semi-analytical and numerical modes. Modes displayed are for 
Q/Q/f = 0.300 (upper line), and Q/Q^ = 0.707 (lower line). Quantum numbers (n,£,m) are as follows: upper line (50,0,0), 
(50, 1,0), (50,0, 1), (50, 1, 1). Lower line (50,0,0), (50, 1,0), (50,2,0), (50,3,0). Blue/black continuous line: numerical results; 
red/dark gray dashed line: semi-analytical results. 



0.15 



0.10 



0.05 



0.00: 




Fig. 6. (Colour online) Frequency spacings S m (€ — 0) = 
<^n,e=o,m ~ w„,f=o,m=o, normalized by a> p , as a function of mj -\Jn. 
The integers n, { and m are the quantum numbers of island 
modes. The rotation rate is £2/n# = 0.419. Blue/black dots: nu- 
merical results. Red/dark gray dashed line: semi-analytical re- 
sults for the last term in Eq. (l56l l. Numerical modes included 
are for quantum numbers n e {42,44,46,48,50}, t — and 
me [-3,3]. 



v = -\fu)(r - ro), where ro is the radial position of the ray y, 
while the value of T is obtained from the eigenvectors of the 
monodromy matrix M. In Fig. [7j the equatorial cuts of semi- 
analytical and numerical modes are plotted for different rota- 
tional velocities and quantum numbers I and m. The chosen 
modes are representative of the different behaviors observed. 
Discrepancies between semi-analytical and numerical results are 
mainly due to edge effects. This occurs when the transverse ex- 
tent of the mode (which scales as 1 / yfuj) reaches either the po- 
lar axis for small rotations, or the surface near the equator for 
high rotations. Finally, avoided crossings can also contravene an 
accurate prediction for mode amplitudes since the amplitudes of 
crossing modes will be linear combinations of all the modes con- 
tributing to the crossing. Hence, modes undergoing an avoided 
crossing can differ significantly from Eq. ( |39l (cf. third panel in 
Fig.UJ. Overall, there is nevertheless a good agreement between 
the semi-analytical and numerical results for mode spatial distri- 
butions, showing the validity of Eq. (139) . 



5. Phenomenology and observables for 
asteroseismology 

In this section, we show that the asymptotic theory provides a 
simple understanding of the evolution of the island mode spec- 
trum with rotation. Then, the physical content of the potentially 
observable frequency spacings 6„, 5( and 5,„ is discussed. 
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Fig. 8. (Colour online) Normalized frequencies q)/u> p as a func- 
tion of normalized rotation velocity Q/£2^ in the observer's 
frame. Different colors correspond to different values of |m|. 
Blue: m = 0, green: |m| = 1, red: \m\ = 2, cyan: \m\ = 3. Modes 
included are for the quantum numbers listed in Table[T](or equiv- 
alentlyn, e [21,25], l s e [0,3],m s e [-£ s ,£ s ]) 



Figure [8] displays the global evolution of all the numerical 
frequencies considered in the observer's frame, whose island 
mode quantum numbers can be found in Table Q] (or equiva- 
lently n s e [21,25], £ s e [0,3], m s e [-£ s ,£ s ] in spherical 
modes quantum numbers). The first phenomenon that can be 
noticed is a global decrease of frequencies with rotation. This 
effect is simply due to the increasing volume of the star when 
it is spinning rapidly. Besides this global effect, the evolution of 
the spectrum's organization can be inferred from the evolution 
of frequency spacings 6„, 5{ and S m . The spacing 6„ stays almost 
constant from null up to high rotations, its value remaining close 
to half the large frequency separation of the spherical model. 
If a large number of island modes are detected in an observed 
spectrum, 6„ should be easily extracted from the data. By con- 
trast, the rapid evolution of 5( with rotation will strongly modify 
the spectrum's organization. This is shown in Fig. [9] where for 
clarity only a few m = modes have been displayed: the £ = 
0, n e [43, 46] and I = 1 , n e [42, 44] modes (or equivalently the 
(n s J s ) € {(21,1), (21, 2), (21, 3), (22,0), (22, 1), (22, 2), (23,0)} 
modes). Starting from the usual structure at zero rotation involv- 
ing the large and small separations of Tassoul's theory, the spec- 
trum reorganization induced by the decrease of 5c can also be 
viewed as an increase of the small separation 5. Then, above 
Q./Q.K — 0.45, the structure of the m = spectrum remains prac- 
tically unchanged. 

Now, to illustrate the evolution of the spectra of different m, 
Fig. [TOl displays the n = 44, t = 0, m e {-2,-1,1,2} mode 
frequencies as a function of the rotation rate together with the 
n e [43,46], I = 0, m = frequencies. The main feature of 
this evolution is the decrease of 6 m from ^\m\ at zero rota- 
tion to very small values at high rotations. When multiple ts of 
island modes are defined as in Sect. [4] and Fig. [3] they show no 
regularity at small rotation rates. Note, however, that the splitting 
o) m —0)- m is always very close to -2mQ because the effects of the 
Coriolis force are negligible. By contrast, at high rotation rates, 
as 6 m vanishes above Q./Q.K — 0.45, the m € {-2,-1,0,1,2} 
modes clearly form a regular multiplet, as can be seen in Fig. [TUJ 
where deviations from strict O spacings are due to the (m/ -\Jn) 2 




0.2 0.3 
Q./Q.K 

Fig. 9. (Colour online) Normalized frequencies a>/co p as a func- 
tion of normalized rotation velocity Q./Q.K- The quantum num- 
bers of the displayed modes are: I — 0, n e [43,46] and £ = 1, 
n e [42,44] with m = 0. Blue: £ = 0, green: I = 1. For clarity, 
the corresponding degrees of the spherical harmonics are also 
written. We have outlined the large frequency separation A, the 
small frequency separation 5 - bi n ,,(\ - <^n,-i/,+2 an d tne <^n> &t 
spacings with arrows. 




0.2 0.3 
Q./Q.K 

Fig. 10. (Colour online) Normalized frequencies cj/oj p as a func- 
tion of normalized rotation velocity £1/Qk- The quantum num- 
bers of the displayed modes are: n = 44, i — 0, m e {-2, -1,1,2} 
with the n e [43,46], £ = 0, m = modes. Blue: m = 0, red: 
\m\ € {1,2}. For clarity, the corresponding degrees of the spher- 
ical harmonics are also written. Continuous line: frequencies in 
the observer frame, dashed line: in the inertial frame. We have 
outlined the large frequency separation A and the 6„, spacing 
with arrows. 



term. Since for such rotation rates the structure of the m = 
spectrum remains unchanged, the evolution of the whole spec- 
trum in the observer's frame is dominated by the advection term 
mfl. 

The global evolution of mode frequencies in the observer's 
frame shown in Fig. [8] also presents some particular events: a 
first clustering of mode frequencies occurs around Q./Q.K — 0.25 
and then a second one around Q./Q.K 0.56. Both phenomena 



10 



M. Pasek et al.: Regular Oscillation Sub-spectrum of Rapidly Rotating Stars 



can be understood from the asymptotic theory. According to the 
asymptotic formulas Eqs. (1411144b . crossings of mode frequen- 
cies will happen when 6c/6„, or equivalently a/71, has a ratio- 
nal value. Though the asymptotic theory predicts true eigenvalue 
crossings, it is known that these crossings wi ll be avoided if the 
two m odes are of the same symmetry class (lLandau & Li fshitz 
119771) . As can be seen on Fig. [4] 5 c becomes equal to 6„ at some 
rotation rate around Q/£2^ 0.25 where the spectrum for a 
given m simplifies to io njn = 6 n (m)n + /3. The degeneracy occurs 
between modes of a different symmetry class, and the rotation 
rate at which it occurs depends only weakly on the m values of 
the modes, if m is small. This property translates itself into a 
clustering of the full spectrum in the observer's frame because 
it turns out that this rotation rate is close to S„/2; and S m , that 
decreases from an initial value of mS n to zero at high rotation, 
is around md n /2 at this intermediate rotation. The second fre- 
quency clustering close to CI/CIk 0.56 is related to the fact 
that d m vanishes at high rotation. In this regime, the different m 
spectra are expected to collapse onto a single spectrum in the ro- 
tating frame but not in the observer's frame. However, when Q 
is equal to 5„, the near degeneracy of the m spectra produces the 
frequency clustering observed at Q./Q.K - 0.56. 

One of the interests of asymptotic theories in asteroseismol- 
ogy is to gain physical insights into seismic observables such as 
6 n , 61, and 5 m . In the following, we briefly discuss this point with 
emphasis on the differences and similarities with the physical 
content of the large and small separation from Tassoul's the- 
ory. The spacing 5„ depends only on the acoustic travel time 
T y — y along the acoustic ray y. We expect T y to be dom- 
inated by the time spent in the sub-surface region where the 
sound speed is much smaller than in the interior. While the path 
of the ray varies with rotation, 6„ remains approximately pro - 
portional to the mean density as shown by Ree se et all ([2008). 
On the other hand, the 6 c spacing depends also on the second 
derivatives of the sound speed transverse to the ray y, an infor- 
mation integrated all along the ray. As long as the path of the 
stable ray goes through the central region of the star, the island 
mode frequencies should be sensitive to the chemical stratifica- 
tion and thus the age of the star. However, after the bifurcation 
at Q./Q.K - 0.26, the ray path progressively avoids the central 
region and the island modes do not contain this information any- 
more. Another interesting property of 6( (or 6c/6 n ) is that it is 
very sensitive to rotation as long as Q/Ojf < 0.35. Finally, for 
high rotation rates (Q./Q.K > 0.40), the value of S m , that can 
be detected through the irregularity of multiplets, also gives an 
information on rotation since it is proportional to tf> ^ds (Eq. 
d56l >) where the distance of the ray to the rotation axis d strongly 
depends on the rotation rate. 

6. Conclusions 

In this paper, we derived an asymptotic formula for frequencies 
that predicts and describes regular spacings in the p-mode spec- 
trum of rapidly rotating stars. The derivation relied on finding 
a stable periodic solution of the acoustic ray dynamics, and ob- 
taining an expression for the modes that are localized around 
this ray, the so-called island modes. The method thus provides a 
formula for the island modes frequencies, as well as a formula 
for the mode spatial distributions. We compared these semi- 
analytical formulas with results from numerical computations of 
high-frequency oscillations in rotating polytropic stellar models. 
The frequency spectrum is characterized by the three spacings 
6„, 6( and 6 m . The agreement was shown to be good for 8„ and 



S( at almost all rotation rates, while 6 m shows significant discrep- 
ancies at low rotation rates. The spacing S n stays almost constant 
at all rotation rates with a value that is close to half the large fre- 
quency separation of the non-rotating model. On the other hand, 
the rapid decrease of S{ strongly modifies the spectrum's orga- 
nization up to Q/f2# =i 0.4, while above that rotation rate, 6c 
remains approximately constant. For such high rotation rates, 
the spacing 8 m nearly vanishes, thus in the observer's frame the 
evolution of the whole spectrum is dominated by the advection 
term m£L. We have also seen that the combined evolution of these 
frequency spacings with rotation leads to particular events such 
as true or near degeneracies, that can significantly simplify the 
spectrum. In addition to these new insights on the evolution of 
the island mode spectrum with rotation, the asymptotic theory 
provides semi-analytical formulas for the regular spacings, in 
particular simple formulas for S n and 6c- 

The present asymptotic theory should be useful for differ- 
ent aspects of stellar seismology in the presence of rapid stellar 
rotation. The regular frequency spacings are potentially observ- 
able and our results provide guidance to look for them in data. 
While investig ations dedicated to the search for regularities are 
necessary (e.g. Lignieres et af] (120101) 1. we expect that the easi- 
est quantities to detect in an island mode spectrum are 8„ at any 
rotation rate, 2mO at small rotation rates, and Q at high rotation. 
For modeling pulsations, the asymptotic theory provides a new 
approach, complementary to numerical computations. One of its 
advantages is to give a quick estimate of frequency spacings for 
a given stellar model, which in turn can be used to search for 
patterns in numerically computed spectra. In the same spirit, the 
semi-analytical amplitude distributions might provide a useful 
approximation for calculating mode visibilities and spectral sig- 
natures. 

The asymptotic theory in itself can be improved and ex- 
tended in various ways. We have already mentioned that the 
method needs to be refined at the rotation rate where the bi- 
furcation of the stable ray occurs, using uniform approxima- 
tions of the ray dynamics. It would also be interesting to predict 
analytically the rotation rate of the bifurcation for a given se- 
quence of stellar models. The present method also assumes that 
the modes are governed by local dynamics around the stable ray. 
This assumption can be tested with a numerical EBK m ethod 
applied to the tori of the stable island (Bohigas et alJll993h . al- 
though this method is complicated to implement in pratice. In 
this paper, we left aside the determination of the actual num- 
ber of modes that are described by the asymptotic theory. An 
estimate of such a number can be obtained by computing sys- 
tematically the p hase space volume of stable islands for differ- 
ent rotations (e.g. Lignie res & Geo rgeot (2009)). Then, knowing 
the value of 6c that gives the mean distance between island mode 
frequencies, or the mean density of these modes, one could com- 
pute the ^ max of the modes that satisfy our formulas. Another 
aspect that we have not modeled is avoided crossing in spite of 
the fact that it will induce important deviations, especially at low 
frequencies. Strong gradients of the sound speed will also pro- 
duce deviations from the asymptotic theory. A technique called 
ray-splittin g, that has already b een used successfully in quan- 
tum chaos (Bliime ret alJll996h . could account for this effect. 
Finally, a similar technique could be applied to asymptotic grav- 
ity modes, that were shown recently to have connections with 
ray theory dBallot et al.ll201ll) . 
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Appendix A: Wave equation in the meridional plane 

In this section, we derive the two-dimensional wave equation 
in the meridional plane of the star. We start from the three- 
dimensional wave equation Eq. <j4j in spherical coordinates 
(r,9,<p), and from the expression of the mode amplitude as 
*P = *P m exp(/ra</>) we obtain the following wave equation: 



A l i d 

K,e + - hr 

r \or 



1 d 



rtan 6 88) (rsin#) 2 



2 2 
U) — U)7, 



Hi j Ivu _ o 

A m ^ -i L m — w , 

C 

(A.l) 

where 9 is the colatitude. We want to cancel out of this equa- 
tion the terms multiplied by first order derivatives of T*,,, in order 
to obtain a two-dimensional Helmholtz-like equation. We thus 
introduce the ansatz *P m (r, f?) = y3(r, 8)<& m (r, 9) in Eq. ( lA.lt to 
obtain that B(r, 9) must satisfy 



Bdr r 



to cancel out terms in and 



dr 



2 dp 



1 



--fr + 



B 88 tan 9 



= 



to cancel out terms in This leads to the solution 



B(r,6) 



B 



(A.2) 



(A.3) 



(A.4) 



Vr sin# 

and we choose B = 1 to yield Eq. (0. 

Appendix B: Expression of the imaginary part of T 

In this section, we derive the expression for the imaginary part 
of r using the Wronskian of Eqs. d23ll24l i. For these equations, 
the Wronskian is defined as 



W(s) = z(s) 



1 dz(s) 1 dz(s) 



c s (s) ds c s {s) ds 



z{s). 



corresponding to 



W(s) = z(s)p(s) - p(s)z(s) 



(B.l) 



(B.2) 



where (z, p) and (z, p) are the two independent solutions of Eqs. 
( 123H24I I. Now, in order to obtain the variation of the Wronskian 
with s, we need the Abel's differential equation identity that says 
that for an equation of the following type 



d 2 v dy 

-ji+P(x)-f + Q(x)y = 0, 
dx L dx 



with the Wronskian 



W: 



dy 2 dyi 
■yi- —yi ■ 



(B.3) 



(B.4) 



dx dx ' 

and y\,y2 the two independent solutions of Eq. ( IB. 31 ). there is the 
identity 

W(x) = Wo exp I - | P(x')dx') . (B.5) 



If we apply this identity to the equation on z(s) that can be de- 
rived from Eqs. (1231124b we obtain that 



W(s) = W(s ) exp 



("I 



Oxc, ds' 



This yields that the Wronskian stays constant for all s: 

W(s) = W(s Q ) = Wo . 
From Eqs. (|2B, dB.lt and ( IB.7t we thus obtain that 



zz 



yielding 



and therefore 



Im[F] = - 



ImrT] = - 



Wn 



ZZ 



Wn 



(B.6) 



(B.7) 



(B.8) 



(B.9) 



(B.10) 



Hence this shows that Im[F] keeps a constant sign along the ray 
when s varies. Thus if T is chosen such that it localizes the func- 
tion at so, then the function will stay localized around the ray for 
all s. Since the localization is obtained for Im[T] > or equiva- 
lently iWo > 0, we choose the Wronskian such that 



W(s) = -i . 



(B.ll) 



Appendix C: Analogy between the function z and 
the deviation of nearby rays 

In this section, we show that the equation satisfied by the func- 
tion z corresponding to Eqs. (!23H24t is the same as the equation 
describing the deviation of two nearby rays. It will thus be pos- 
sible to compute the evolution of z from the ray dynamics. To be 
an acoustic ray of the system, the ray y must correspond to an 
extremum of the action S defined as 

rif rif \ 

S(qi,qf,co) = I pdq= —dcr. (C.l) 

Jqi Jq, c s 

The length element dcr can be expressed in the coordinate sys- 
tem centered on the ray y defined in Sect. 13. li as 



dcr 2 = (1 -£k(s)) 2 ds 2 +df , 



yielding 



do- 2 = 



ds 2 



(C2) 



(C3) 



The action function S can therefore be expressed in this coordi- 
nate system as 



ds. 



(C4) 



with the corresponding Lagrangian 
1 



(W*(*)) 2 + |^ 



(C5) 
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The Lagrangian L becomes for £ small: 



1 



Us) 



1 



1 dc s 



2 c/i) <9£ 



k(s) dc s 2 1 d 2 c 



i >2 



g,(j) d% 2c s (s) df 



(C6) 



1 l<*tfg 



+ 0(f) 



where £ = d^/ds and derivatives in £ are evaluated at £ = 0. As is 
known from classical mechanics, for the action to be extremum, 
the associated Lagrangian must satisfy the Euler-Lagrange equa- 
tion 

Keeping only quadratic terms in £ from this equation yields 



c s (s) afs 

i 



1 



c s (s) 2 \ a? / ' c s (s) d% 



2c s (s) de 



«(*)• 



1 <9c s 



(C8) 



+ 0(^ 3 ) = . 



Us) d$ 

This equality must be valid for all £ so we can deduce that: 

1 dc s 



K(S) 



Us) d% 



(C9) 



and 



1 



1 d% 



c s (s) ds\c s (s) dsj Us) 2 



1 IdU 2 



+ 



k(s) dc s 

U.s)- \ <)fl ' Us~)~di 



2~us) oe 



€ = 0. 

(CIO) 



We thus obtain the equation for the deviation between nearby 
rays as 

1 d I 1 d%\ 



U s ) ds \ U S ) ds 



+ K(s)% = . 



(Cll) 



and this equation is the same as the equation satisfied by z that 
one can derive from Eqs. (1231124b . 

Appendix D: Eigenvectors and eigenvalues of the 
monodromy matrix M 

In this section, we derive the necessary formulas to express the 
stability angle a and the function T in terms of the elements of 
the monodromy matrix M. The characteristic polynomial of a 
2x2 matrix M is 



A - Tr( [M)A + det(M) = . 



(D.l) 



Then, it is known that a monodromy matrix is a symplectic ma- 
trix, which implies that det(M) = 1. For a stable ray, we have 
|Tr(M)| < 2, and thus the roots of this polynomial become 



Tr(M) y/-Tr(M) 2 + 4 
*U = — ±» 2 ■ 

We thus obtain the expression for the stability angle a as 



(D.2) 



exp(+/a) = - [Tr(M) ± i V-Tr(M) 2 + 4] , (D.3) 



and we obtain Eq. ( l34l >. 

The eigenvector Vi will then be given by the following equa- 
tion: 

(Mn-Ai Mi 2 \(zi s 



M 2 i M 22 -Ai]\p\ 
which corresponds to 



= 0, 



(M n - Ai)z\ +M x2 p\ = 
M 2 \Zi + (M 22 - Ai)px = . 

Recalling that F(s) = p(s)/z(s) we write 

Zl 

and thus we obtain that 

(Mn - AO M 2i 



M u (M 22 - Ai) 



(D.4) 



(D.5) 
(D.6) 



(D.7) 



(D.8) 



This yields 



_ -M u + M 22 . V-Tr(M) 2 +4 

1 — ^ + I ^ , \\J.y) 



2M 12 2Mn 
and with the same procedure we obtain the expression of T2: 

-M n +M 22 . V-Tr(M) 2 + 4 



2M 
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